Development of an Assessment Model for the Effect of the Replacement of Minimal Artificial Ossicles on Hearing in the Inner Ear

Due to ethical issues and the nature of the ear, it is difficult to directly perform experimental measurements on living body elements of the human ear. Therefore, a numerical model has been developed to effectively assess the effect of the replacement of artificial ossicles on hearing in the inner ear. A healthy volunteer’s right ear was scanned to obtain CT data, which were digitalized through the use of a self-compiling program and coalescent Patran-Nastran software to establish a 3D numerical model of the whole ear, and a frequency response of a healthy human ear was analyzed. The vibration characteristics of the basilar membrane (BM) after total ossicular replacement prosthesis (TORP) implantation were then analyzed. The results show that although the sound conduction function of the middle ear was restored after replacement of the TORP, the sensory sound function of the inner ear was affected. In the low frequency and medium frequency range, hearing loss was 5.2~10.7%. Meanwhile, in the middle–high frequency range, the replacement of a middle ear TORP in response to high sound pressure produced a high acoustic stimulation effect in the inner ear, making the inner ear structures susceptible to fatigue and more prone to fatigue damage compared to the structures in healthy individuals. This developed model is able to assess the effects of surgical operation on the entire hearing system.


Introduction
Since the beginning of the 21st century, the two journals "Science" and "Nature" have respectively reported that "The problems of deafness are deeper and more complex, which are even more important than those of blindness" and "Restoration of auditory function remains a lofty, (challenge), but achievable goal for deaf patients and for scientists. Because hearing loss is a multifactorial problem arising from many possible causes, it is unlikely that any single strategy will be an effective cure-all" [1]. Thus, up to now, the investigation of hearing mechanisms in humans is still a medical challenge for mankind [2].
Middle ear lesions can lead to conductive deafness. At present, artificial ossicular replacement is an effective method for the treatment of conductive deafness. Some research has been conducted in this relevant area. Murphy [3] investigated and evaluated the degree of hearing recovery after total ossicular replacement prosthesis (TORP) implantation. Yao et al. studied the effect of artificial ossicular replacements with different materials on postoperative hearing recovery [4]. Kelly et al. [5] investigated the effect of different artificial ossicular quality on sound transmission. Marchese  The geometric dimensions of the numerical model were based on images obtained from CT scanning of healthy human ears at Zhongshan Hospital (GE lightspeed VCT 64 slice spiral CT machine; scanning parameters: collimation 0.625 mm, tube rotation time 0.4 s, reconstruction layer thickness 0.625 mm, interval 0.5-0.625 mm, and a voxel size of 1.0 × 1.0 × 0.75 mm 3 ). CT scanning images were digitized by a self-compiling program and then imported into the Patran software to reconstruct the three-dimensional model of the human ear. The mesh was then divided and the boundary conditions and material parameters were set. Finally, three-dimensional fluid-solid coupling dynamic analysis was carried out using the Nastran software. This study was approved by the Ethics Committee of the Zhongshan Hospital affiliated to Fudan University. Written informed consent was obtained from the patient.
The geometric dimensions of the numerical model were based on images obtained from CT scanning of healthy human ears at Zhongshan Hospital (GE lightspeed VCT 64 slice spiral CT machine; scanning parameters: collimation 0.625 mm, tube rotation time 0.4 s, reconstruction layer thickness 0.625 mm, interval 0.5-0.625 mm, and a voxel size of 1.0 × 1.0 × 0.75 mm 3 ). CT scanning images were digitized by a self-compiling program and then imported into the Patran software to reconstruct the three-dimensional model of the human ear. The mesh was then divided and the boundary conditions and material parameters were set. Finally, three-dimensional fluid-solid coupling dynamic analysis was carried out using the Nastran software. This study was approved by the Ethics Committee of the Zhongshan Hospital affiliated to Fudan University. Written informed consent was obtained from the patient.

Material Properties
The material properties of the ossicular chain numerical model are shown in Table 1 [19][20][21][22][23][24], and the material properties of the soft tissue finite element model (FEM) are shown in Table 2 [24,25]. The Poisson's ratio of each part of the middle ear structure is 0.3, the structural damping coefficient is 0.4, the viscosity of the fluid is 0.001 NS/m 2 , and the damping coefficient β of the fluid is 0.0001 s.

Material Properties
The material properties of the ossicular chain numerical model are shown in Table 1 [19][20][21][22][23][24], and the material properties of the soft tissue finite element model (FEM) are shown in Table 2 [24,25]. The Poisson's ratio of each part of the middle ear structure is 0.3, the structural damping coefficient is 0.4, the viscosity of the fluid is 0.001 NS/m 2 , and the damping coefficient β of the fluid is 0.0001 s.  The material properties of the inner ear structure shown above were obtained from the relevant published references [26][27][28]. The material properties of each part of the inner ear in the numerical model in this paper are as follows: Oval window: the elastic modulus is E = 0.2 MPa, Poisson's ratio is µ = 0.3, density is ρ = 1200 kg/m 3 , and the damping coefficient is β = 0.5 × 10 −4 s. Round window: the elastic modulus is E = 0.35 MPa, Poisson's ratio is µ = 0.3, and the damping coefficient is β = 0.5 × 10 −4 s. Lymphatic fluid (scala vestibuli, scala tympani, scala media, 3 semicircular canals, and lymphatic fluid in the vestibuli): density is ρ = 1000 kg/m 3 , sound velocity is C = 1400 m/s, the damping coefficient is β = 1.0 × 10 −4 s, and viscous damping is D = 0.001 NS/m. BM: As the length of the BM changes, the elastic modulus decreases linearly from 50 MPa at the base of the cochlea to 15 MPa at the middle and then decreases linearly to 3 MPa at the apex. The damping coefficient β varies linearly from 0.2 × 10 −3 s at the base to 0.1 × 10 −2 s at the apex, with a Poisson's ratio of 0.3.

Boundary Conditions of the Whole Ear FEM
Because the boundary conditions need to be provided for the numerical model's calculations according to the structural characteristics of the human ear and the connection relationship and related characteristics between the structures of the human ear, the boundaries of some tissues of the human ear are set appropriately based on the mechanical principle. The details are as follows: (1) Application of 80 dB SPL (0.2 Pa), 90 dB SPL (0.632 Pa), and 105 dB SPL (3.56 Pa) surface pressure to the opening surface of the external ear canal or TM to simulate pure tone sound pressure stimulation (100-10,000 Hz); (2) The positions of soft tissues (tensor tympani, superior mallear ligaments, anterior mallear ligaments, lateral mallear ligament, superior incudal ligament, posterior incudal ligament, stapedial tendon) associated with the temporal bone were defined as the fixed constraint (constrain all displacement and all rotation); (3) The outer edge of the TM's annular ligament was defined as the hinged constraint (only constrains all displacement, not rotation); (4) The outer edge of the SF annular ligament was defined as the fixed constraint (constrain all displacement and all rotation); (5) The outer edge of the oval window and the round window were fixed constraints; (6) The three edges of the BM (both sides and the base of the cochlea) were considered as hinged constraints (only constrains all displacement, not rotation); (7) The external ear canal wall and the inner ear bony labyrinth wall were set as the rigid wall;

Establishment of the TORP Model
In this study, an artificial ossicular TORP (produced by the Medtroni pany in the USA) was applied for FE analysis. The size of the TORP (363) w the length was 6.7 mm in total and the diameter of its columella was 0.9 mm the undersurface diameter of the round terminal disc was 3.2 mm, and the e columella was a plane. The model and mesh of the TORP were established tran FE software. The mesh was divided into 539 nodes (100 eight-node hexa elements and 325 six-node pentahedral (Wedge6) elements) and the elemen defined as Solid. The round terminal disc of the TORP was connected to th TM and its columella's end face was connected to the SF or oval window m ditionally, a 0.5 mm cartilage slice was set between the TM and the under round terminal disc of the TORP. The material properties of the titanium were as follows [5]: the elastic modulus was 116 GPa, density was 4500 k Poisson's ratio was 0.33. The mesh division of the model is shown in Figure   Figure 3. Whole ear FEM after TORP replacement.

Establishment of the TORP Model
In this study, an artificial ossicular TORP (produced by the Medtronic Xomd Company in the USA) was applied for FE analysis. The size of the TORP (363) was as follows: the length was 6.7 mm in total and the diameter of its columella was 0.9 mm. In addition, the undersurface diameter of the round terminal disc was 3.2 mm, and the end face of the columella was a plane. The model and mesh of the TORP were established using the Patran FE software. The mesh was divided into 539 nodes (100 eight-node hexahedral (Hex8) elements and 325 six-node pentahedral (Wedge6) elements) and the element property was defined as Solid. The round terminal disc of the TORP was connected to the center of the TM and its columella's end face was connected to the SF or oval window membrane. Additionally, a 0.5 mm cartilage slice was set between the TM and the undersurface of the round terminal disc of the TORP. The material properties of the titanium in the TORP were as follows [5]: the elastic modulus was 116 GPa, density was 4500 kg·m −3 , and the Poisson's ratio was 0.33. The mesh division of the model is shown in Figure 3.

Establishment of the TORP Model
In this study, an artificial ossicular TORP (produced by the pany in the USA) was applied for FE analysis. The size of the TOR the length was 6.7 mm in total and the diameter of its columella w the undersurface diameter of the round terminal disc was 3.2 mm columella was a plane. The model and mesh of the TORP were e tran FE software. The mesh was divided into 539 nodes (100 eightelements and 325 six-node pentahedral (Wedge6) elements) and th defined as Solid. The round terminal disc of the TORP was conne TM and its columella's end face was connected to the SF or oval w ditionally, a 0.5 mm cartilage slice was set between the TM and round terminal disc of the TORP. The material properties of th were as follows [5]: the elastic modulus was 116 GPa, density w Poisson's ratio was 0.33. The mesh division of the model is shown   The frequency response curves of the umbo and SF were obtained through numerical simulation analysis and compared to the experimental data of Gan et al. [29,30], as shown in Figures 4 and 5. external ear canal.

Load Acting on the TM
The frequency response curves of the umbo and SF were ob simulation analysis and compared to the experimental data of G in Figure 4 and 5.   The frequency response curves of the umbo and SF were ob simulation analysis and compared to the experimental data of G in Figure 4 and 5.

Verification of the FEM of the Spiral BM
Since the vibration of external sound through the TM and the ossicular chain drives the vibration of the SF, it is the vibration of the SF that provokes the response of the BM through cochlear lymphatic fluid, and the vibration of the stapes has a direct relationship with the response of the BM. The ratio of the amplitude of the BM to the amplitude of the SF represents the response of the BM caused by the vibration of the SF per unit, which clearly reflects the amplification characteristics and frequency selection characteristics of the BM.
Considering the influence of the spiral BM in the cochlea, a sound pressure of 90 dB SPL was applied to the opening surface of the external ear canal for frequency response analysis, as shown in Figure 2.
The BM data derived from the spiral inner ear numerically simulated by this model, and the experimental data at 12 mm from the base of the BM obtained by Gundersen et al. and Stenfelt et al. [12,13], were compared to some extent, and the BM to SF amplitude ratio-frequency response curves were constructed, as shown in Figure 6. chines 2023, 14, x FOR PEER REVIEW SF represents the response of the BM caused by the vibration of the SF clearly reflects the amplification characteristics and frequency selection c the BM.
Considering the influence of the spiral BM in the cochlea, a sound pr SPL was applied to the opening surface of the external ear canal for freq analysis, as shown in Figure 2.
The BM data derived from the spiral inner ear numerically simulated and the experimental data at 12 mm from the base of the BM obtained b al. and Stenfelt et al. [12,13], were compared to some extent, and the BM t ratio-frequency response curves were constructed, as shown in Figure 6.  As shown in Figure 6, in the frequency range of 100-1000 Hz, the amplitude ratio-frequency response curve of the BM and SF obtained by numerical simulation with this model is slightly higher than that of Gundersen  In the frequency range of 3000-8000 Hz, the amplitude ratio-frequency response curve of the BM to the SF obtained by numerical simulation is close to the experimental curve of Gundersen et al., and the curve trend is consistent with the peak value of this frequency range. The peak value appears at about 3500 Hz, and the average error in this interval is about 6%.
In Figure 6, comparing these curves, the amplitude ratio-frequency response curves of the BM to the SF measured by Gundersen et al. and Stenfelt et al. have two peaks, and those derived from numerical simulations of the spiral cochlea in the model also have two peaks over the entire frequency range (100-8000 Hz), one at approximately 400 Hz and the other at 3500 Hz, which occur in the low and mid-frequency ranges, respectively. This shows that the vibration characteristics of the BM obtained from FE method simulation of the spiral inner ear (including the spiral BM) are closer to those obtained from inner ear experiments, that is, FE simulation of the spiral inner ear may also accurately reflect the inner ear characteristics of the actual human inner ear, thus providing good theoretical guidance and help for clinical research in this area. Figure 7 more clearly shows the error in the results of the model calculation in this paper and experimental results in terms of frequency. In the frequency range of 1000-3000 Hz, the average error was minimal (only 1.6%). In the frequency range of 100-1000 Hz, the average error was 12%. In the frequency range of 1000-10,000 Hz, the average error was 6%. The overall average error was 8%.
achines 2023, 14, x FOR PEER REVIEW 1.6%). In the frequency range of 100-1000 Hz, the average error was 12%. In the f range of 1000-10,000 Hz, the average error was 6%. The overall average error w In summary, the SF and umbo data obtained from the human ear FEM sim in this paper are close to the experimental data in terms of amplitude and trend tively, thus verifying the correctness of this model. In addition, the simulation vibrations on the spiral BM and a comparison with experimental data were used the correctness of the spiral cochlea (contains the spiral BM).

Numerical Analysis of the Vibration Characteristics of the BM after Replacement w TORP
In this paper, the frequency response analysis of the FEM of the whole hea tem after replacement with a titanium TORP in the frequency range of 100-10,00 performed by applying a sound pressure of 90 dB SPL over the opening surfa external ear canal while considering the effect of the spiral BM built into the co shown in Figure 8. Since the vibration of external sound through the TM and the ossicular cha the vibration of the SF, it is the vibration of the SF that provokes the response o through cochlear lymphatic fluid, and the vibration of the stapes has a direct rel with the response of the BM. In order to clearly reflect the amplification chara and frequency selection characteristics of the spiral BM after replacement with a In summary, the SF and umbo data obtained from the human ear FEM simulations in this paper are close to the experimental data in terms of amplitude and trend, respectively, thus verifying the correctness of this model. In addition, the simulation of nodal vibrations on the spiral BM and a comparison with experimental data were used to verify the correctness of the spiral cochlea (contains the spiral BM).

Numerical Analysis of the Vibration Characteristics of the BM after Replacement with a TORP
In this paper, the frequency response analysis of the FEM of the whole hearing system after replacement with a titanium TORP in the frequency range of 100-10,000 Hz was performed by applying a sound pressure of 90 dB SPL over the opening surface of the external ear canal while considering the effect of the spiral BM built into the cochlea, as shown in Figure 8.
Micromachines 2023, 14, x FOR PEER REVIEW 1.6%). In the frequency range of 100-1000 Hz, the average error was 12%. In range of 1000-10,000 Hz, the average error was 6%. The overall average er In summary, the SF and umbo data obtained from the human ear FE in this paper are close to the experimental data in terms of amplitude and tively, thus verifying the correctness of this model. In addition, the simul vibrations on the spiral BM and a comparison with experimental data were the correctness of the spiral cochlea (contains the spiral BM).

Numerical Analysis of the Vibration Characteristics of the BM after Replacem TORP
In this paper, the frequency response analysis of the FEM of the who tem after replacement with a titanium TORP in the frequency range of 100performed by applying a sound pressure of 90 dB SPL over the opening external ear canal while considering the effect of the spiral BM built into shown in Figure 8.  Since the vibration of external sound through the TM and the ossicular chain drives the vibration of the SF, it is the vibration of the SF that provokes the response of the BM through cochlear lymphatic fluid, and the vibration of the stapes has a direct relationship with the response of the BM. In order to clearly reflect the amplification characteristics and frequency selection characteristics of the spiral BM after replacement with a titanium TORP, the ratio of BM amplitude to the amplitude of the SF is used in this paper, which indicates the BM response caused by the vibration of the SF.
The amplitude-frequency response curves of the nodes at the center of the SF and at 12 mm from the base of the cochlea on the BM were obtained by simulation analysis, and the relationship between the ratio of the BM amplitude and the stapes amplitude at 12 mm from the base of the cochlea on the BM and frequency was calculated. The BM data were also compared to those of the healthy human ear, as shown in Figure 9. From Figure 9, it can be seen that, in the frequency range of 20-10,000 Hz, th tude ratio-frequency response curve of the BM and the SF after replacement with is closer to the vibration characteristics of the BM of the healthy ear in terms of am and overall trend. There are also two peaks at 600 Hz and 4000 Hz, and the pea amplitude ratio-frequency response curve of the BM to the SF after replacemen TORP is somewhat more lagged than in healthy persons, where the two peaks approximately 400 Hz and 3500 Hz.
In the frequency range of 100-600 Hz, the amplitude ratio of the BM to th creased gradually with frequency after the TORP replacement, and the curve was lower than that of healthy persons. In this frequency range, the difference betw amplitude ratio of the BM and the SF was about 1.03 dB to 3.72 dB, with an avera tive error of about 5.2% compared with that of healthy persons.
In the frequency range of 600-2000 Hz, the BM to SF amplitude ratio after T placement tended to decrease slowly with frequency, and the curve was somewha than that of healthy persons. In this frequency range, the difference between the SF amplitude ratios compared to healthy persons was about 0.67 dB to 5.33 dB, average relative error of about 8.3%.
In the frequency range of 2000-4000 Hz, the BM to SF amplitude ratio afte replacement tended to increase slowly with frequency, and the curve was so lower than that of healthy persons. In this frequency range, the difference betw amplitude ratio of the BM and the SF was about 2.01 dB to 5.54 dB compared wit From Figure 9, it can be seen that, in the frequency range of 20-10,000 Hz, the amplitude ratio-frequency response curve of the BM and the SF after replacement with a TORP is closer to the vibration characteristics of the BM of the healthy ear in terms of amplitude and overall trend. There are also two peaks at 600 Hz and 4000 Hz, and the peak of the amplitude ratio-frequency response curve of the BM to the SF after replacement with a TORP is somewhat more lagged than in healthy persons, where the two peaks occur at approximately 400 Hz and 3500 Hz.
In the frequency range of 100-600 Hz, the amplitude ratio of the BM to the SF increased gradually with frequency after the TORP replacement, and the curve was slightly lower than that of healthy persons. In this frequency range, the difference between the amplitude ratio of the BM and the SF was about 1.03 dB to 3.72 dB, with an average relative error of about 5.2% compared with that of healthy persons.
In the frequency range of 600-2000 Hz, the BM to SF amplitude ratio after TORP replacement tended to decrease slowly with frequency, and the curve was somewhat higher than that of healthy persons. In this frequency range, the difference between the BM and SF amplitude ratios compared to healthy persons was about 0.67 dB to 5.33 dB, with an average relative error of about 8.3%.
In the frequency range of 2000-4000 Hz, the BM to SF amplitude ratio after TORP replacement tended to increase slowly with frequency, and the curve was somewhat lower than that of healthy persons. In this frequency range, the difference between the amplitude ratio of the BM and the SF was about 2.01 dB to 5.54 dB compared with that of healthy persons, with an average relative error of about 10.7%.
In the frequency range of 4000-7000 Hz, the amplitude ratio of the BM to the SF after TORP replacement decreased linearly with frequency, and the curve was higher than that of healthy persons. In this frequency range, compared with the amplitude ratio of the BM to the SF in healthy persons, the difference was about 1.89-9.22 dB, and the average relative error was about 9.1%.
In the frequency range of 7000-10,000 Hz, the amplitude ratio of the BM to the SF after TORP replacement was lower than that of healthy persons. At this frequency range, the amplitude ratio of the BM to the SF was about 1.32-6.43 dB different from that of healthy persons.
A fluid-solid coupling cloud diagram of the whole hearing system (ossicular chain, a displacement diagram of the BM, and a pressure diagram of the external ear canal air and inner ear perilymph) is shown in Figures 10 and 11.  Through numerical simulation with this model, the displacement cloud diagram of the BM after TORP replacement at different frequencies (500 Hz, 1000 Hz, 4000 Hz) is compared with that of healthy human ears, as shown in Figure 12 and 13 below.  Through numerical simulation with this model, the displacement cloud diagram of the BM after TORP replacement at different frequencies (500 Hz, 1000 Hz, 4000 Hz) is compared with that of healthy human ears, as shown in Figure 12 and 13 below. Through numerical simulation with this model, the displacement cloud diagram of the BM after TORP replacement at different frequencies (500 Hz, 1000 Hz, 4000 Hz) is compared with that of healthy human ears, as shown in Figures 12 and 13 below. chain in middle ear; (b) BM's displacement diagram (visual angle 1); (c) BM's displacement diagram (visual angle 2).
Through numerical simulation with this model, the displacement cloud diagram of the BM after TORP replacement at different frequencies (500 Hz, 1000 Hz, 4000 Hz) is compared with that of healthy human ears, as shown in Figure 12 and 13 below.  From the comparison of Figure 12 and 13, the following can be determined: At 500 Hz, the displacement of the BM in the healthy human ear ranged from 0 to 7.63 × 10 −6 mm, and the maximum value of this displacement was observed near the base of the cochlea. The displacement of the BM in the human ear after TORP replacement ranged from 0 to 7.48 × 10 −6 mm, and the maximum displacement value was observed near the apex of the cochlea. BM displacement after TORP was smaller than in the healthy human ear, with the maximum displacement occurring at different positions; however, the values of the displacements were in the same order of magnitude.
At 1000 Hz, the displacement of the BM in the healthy human ear ranged from 0 to 4.83 × 10 −6 mm, and the maximum value of displacement at this time was close to the apex of the cochlea. Displacement of the BM in the human ear after TORP replacement ranged from 0 to 2.65 × 10 −6 mm, and the maximum displacement value was observed close to the apex of the cochlea. The position of the maximum displacement at this frequency is closer to that of the BM in the healthy human ear, though the values are somewhat different. Compared to the healthy human ear, BM displacement after TORP replacement decreased.
At 4000 Hz, the displacement of the BM in the healthy human ear ranged from 0 to From the comparison of Figures 12 and 13, the following can be determined: At 500 Hz, the displacement of the BM in the healthy human ear ranged from 0 to 7.63 × 10 −6 mm, and the maximum value of this displacement was observed near the base of the cochlea. The displacement of the BM in the human ear after TORP replacement ranged from 0 to 7.48 × 10 −6 mm, and the maximum displacement value was observed near the apex of the cochlea. BM displacement after TORP was smaller than in the healthy human ear, with the maximum displacement occurring at different positions; however, the values of the displacements were in the same order of magnitude.
At 1000 Hz, the displacement of the BM in the healthy human ear ranged from 0 to 4.83 × 10 −6 mm, and the maximum value of displacement at this time was close to the apex of the cochlea. Displacement of the BM in the human ear after TORP replacement ranged from 0 to 2.65 × 10 −6 mm, and the maximum displacement value was observed close to the apex of the cochlea. The position of the maximum displacement at this frequency is closer to that of the BM in the healthy human ear, though the values are somewhat different. Compared to the healthy human ear, BM displacement after TORP replacement decreased.
At 4000 Hz, the displacement of the BM in the healthy human ear ranged from 0 to 1.49 × 10 −7 mm, and the maximum displacement value was observed in the first lap of the upper cochlea. The displacement of the BM in the human ear after TORP replacement ranged from 0 to 1.17 × 10 −7 mm, and the position of maximum displacement at this frequency was basically the same as it was in the healthy human ear. BM displacement after TORP replacement was smaller than it was in the healthy human ear, though the values of the displacements were in the same order of magnitude.

Conclusions
After the TORP replacement based on the whole hearing system, the vibration characteristics of the spiral BM changed at low and high frequencies. Compared to the healthy human ear, the position of the peak amplitude value in the BM after TORP replacement was somewhat more lagged than in healthy persons based on the amplitude-frequency response curve. In the low frequency range, BM amplitude after TORP replacement was lower than in the healthy human ear, and the decrease increased as frequency increased, with an average relative difference between the two of 5.2%. In the middle frequency range, BM amplitude after TORP replacement increased slowly with increasing frequency, amplitude was lower than in the healthy human ear, and the average relative difference between the two was 10.7%. In the middle-high frequency range, BM amplitude after TORP replacement is higher than in the healthy human ear, and the average relative difference between the two was 9.1%. In the high frequency range, BM amplitude after TORP replacement was lower than in the healthy human ear, and the average relative difference between the two was 8.6%.
In conclusion, TORP replacement in the middle ear also has some infaust effects on the acoustic sensing function of the inner ear, especially in the low frequency and medium frequency range. Meanwhile, in the middle-high frequency range, TORP replacement in the middle ear in response to high sound pressure produces a high acoustic stimulation effect on the inner ear, making the inner ear structures susceptible to fatigue and more prone to fatigue damage.
The numerical model in this paper can not only study the level of sound conduction recovery in the human middle ear, but also predict the effect of clinical surgery involving TORP replacement on the whole hearing system.